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The potential for neutrino telescopes to discover charged stau production in neutrino-nucleon 
interactions in Earth depends in part on the stau lifetime and range. In some supersymmetric 
scenarios, the next lightest supersymmetric particle is a stau with a decay length on the scale of 
10 km. We evaluate the electromagnetic energy loss as a function of energy and stau mass. The 
energy loss parameter j3 scales as the inverse stau mass for the dominating electromagnetic processes, 
photonuclear and e^e“ pair production. The range can be parameterized as a function of stau mass, 
initial energy and minimum final energy. In comparison to earlier estimates of the stau range, our 
results are as much as a factor of two larger, improving the potential for stau discovery in neutrino 
telescopes. 


I. INTRODUCTION 


Neutrino telescopes have great potential for probing 
physics beyond the Standard Model when new parti¬ 
cles are created by high energy cosmic neutrino passing 
through the Earth. A recent proposal Q showed that the 
detection of two nearly parallel charged particle tracks 
from a pair of metastable staus, which are produced as 
secondary particles in very high energy neutrino-nucleon 
interactions, may provide a unique way of probing the 
supersymmetry (SUSY) breaking scale in weak scale su¬ 
persymmetry models. 

Ultrahigh energy neutrinos from astrophysical sources 
travel undeflected to the Earth. They would then inter¬ 
act with the nucleons in the Earth and pair produce two 
superparticles via the t-channel chargino or neutralino 
exchanges: vN Iq. In SUSY models where grav- 
itino is the lightest supersymmetric particle (ESP) and 
stau is the next lightest supersymmetric particle (NLSP), 
heavier slepton and squarks instantly decay into the stau 
NLSP 12 • The high energy stau would then travel a long 
distance, typically more than tens to 10^ kilometers, be¬ 
fore it decays into the gravitino LSP. The superparticle 
pair production rate, however, is very small in compari¬ 
son to the standard model charged-current cross section. 
Thus production of staus from downward neutrinos is 
negligible compared to showers from production through 
the standard model processes. However, the observability 
of the stau produced via upward going neutrinos could 
be enhanced. The small stau production cross section 
are compensated by the enlarged effective volume of an 
instrumented region due to the very long lifetime of the 
stau and its relatively slow energy loss in the water, ice 
or rock surrounding the detector [J. 

The potential for neutrino telescopes to discover stau 
production in neutrino-nucleon interactions in Earth de¬ 
pends in part on the stau lifetime and range. In order to 
refine estimates of the signal from metastable staus, it is 
crucial to determine the energy loss of the high energy 
stau as it traverses the Earth. The average energy loss 


of a particle traversing a distance X is given by 

-(f)=a + /3E (1) 

where E is the energy of the particle, a describes ion¬ 
ization energy loss, and (3 is the radiative energy loss. 
Ref. Q estimated the stau range using Eq. o with 
the approximation that both a and ff are constants with 
respect to energy. In this paper, we evaluate both the 
energy and mass dependence of the stau electromagnetic 
energy loss. Our results are applicable to any metastable 
heavy charged slepton. We determine f3 and we evaluate 
the average stau range using a one-dimensional Monte 
Carlo simulation of stau propagation in rock. We com¬ 
ment on the implications of the new results for the evalua¬ 
tion of event rates by Albuquerque, Burdman and Chacko 
in Ref. 0. 

This paper is organized as follows. In Sec. nn we an¬ 
alyzed in detail the stau energy loss through photonu¬ 
clear scattering, pair production and bremsstrahlung. 
In Sec. EH we show the stau range using both a one¬ 
dimensional Monte Carlo evaluation and a parametrized 
approximate formula. In Sec. nvi we conclude. 


II. STAU ENERGY LOSS 


The radiative energy loss (3 receives contributions from 
bremsstrahlung, pair production and photonuclear scat¬ 
tering: 


N 

P\E) = - 



dy y 


dcF'’{y,E) 

dy 


( 2 ) 


where y is the fraction of stau energy loss in the radiative 
interaction: 


y = 


E-E' 

E 


( 3 ) 


for final stau energy E'. The superscript i denotes 
bremsstrahlung (brem), pair production (pair) and pho- 
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FIG. 1: Parameter space in (mj-, y/F) in which the stau range 
is dominated by either ionization energy loss or the lifetime, 
when PE a. 



E [GeV]] 


tonuclear (nuc) processes. Avogadro’s number is N and 
the atomic mass number of the target nucleus is A. 

As we show below, the value for [3 decreases with in¬ 
creasing stau mass. The ionization energy loss parameter 
a is nea:^ constant as a function of mass of the particle, 
namely Q 


a ~ 2 X 10"^ GeVcmVg • (4) 

For low enough energies, where [3E <C a (for P ^ 5x 10“® 
cm^/g-(150 GeV/mj), it corresponds to A «C 4 x 10® 
GeV-(mj/150 GeV)), either the lifetime or ionization en¬ 
ergy loss determines the stau range, which scales linearly 
with energy. This is shown in Fig. ^ These two re¬ 
gions are separated by the line which corresponds to 
CTp = mjja, where mj is the stau mass, p is the den¬ 
sity, taken to be the density of standard rock p = 2.65 
g/cm®, and 


cr = 




4 / \ 5 

100 GeV \ 


10^ GeV 


10 km 


( 5 ) 


is the stau lifetime. The quantity \/F is the supersymme¬ 
try breaking scale. In low energy SUSY breaking mod¬ 
els, the SUSY breaking scale y/F could be as low as 10^ 
GeV[2, leading to a stau decay length of 10 km or longer. 


A. Stau Energy Loss due to Photonuclear 
Interactions 

We obtain a photonuclear cross section that describes 
the interactions of charged staus (f, or more generically 
charged sleptons 1) with nuclei via virtual photon ex¬ 
change. Our approach 0 is to treat the IN IX using 
the deep-inelastic scattering formalism in terms of the 
two nucleon electromagnetic structure functions Fi and 

F2. 


FIG. 2: Energy loss of stau due to photonuclear (nuc) interac¬ 
tions, pair production (pair) and bremsstrahlung (brem) for 
a stau mass of 150 GeV. 
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FIG. 3: Energy loss of stau due to photonuclear and pair 
production for rrif = 50 GeV, rtif = 150 GeV and = 250 
GeV. 


We use F 2 parameterized by Abramowitz, Levin, Levy 
and Maor (ALLM) to fit HERA data for a wide range 
of and x. The ALLM parameterization of the struc¬ 
ture function is based on the Regge approach and con¬ 
sists of a Pomeron and a Reggeon contribution 0. Other 
approaches 0, which focus on the diffractive region, but 
also include a perturbative component, give qualitatively 
similar results for the photonuclear energy loss of leptons. 

For stau scattering with nucleons, f{k)N{p) 
T{k')X, the differential cross section is given by 


(Pa{x, Q^) 
dQ^ dx 


47ra^ F2{x, Q^) 

X 


1 - y + 


X 


( 6 ) 
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2/^(1 + /Q'^) 

4(l + i?(a:,Q2)) 


where —Q^ = = (^k — k'Y, x = Q‘^l{2p- q), y as in 

Eq. 0, mj is the stau mass, and M is the nucleon mass. 
R{x, Q^) parametrizes the violation of the Callan-Gross 
relation 2xFi{x,Q'^) = F 2 {x,Q'^): 


R{x, 


{l + ^^)F2{x,Q^)-2xF,ix,Q^) 
2xFi{x, Q2) 


(7) 


Since R{x,Q‘^) is small we approximate R — 0. 

To obtain /3, we convert the differential cross section 
to a dependence on y and Q^: 


ffajy.Q'^) ^ x cPa{x,Q'^) 
dQ'^ dy y dQ'^ dx 

We use the following limits of integration: 

QLn < Q^<2MEy-{{M + m^)^-M^) (9) 

2/min < y<l-mi/E, (10) 

where ~ rnfy'^l{\ - y) and j/min - {{M + - 

M^)/{2ME). 

For comparison, the photonuclear cross section for the 
standard model partner, the charged lepton, is given by 

i 


d^a{x, Q^) 
dQ^ dx 


E 2 {x,Q'^) Mxy 

~ ~y~ 2E 


( 11 ) 


/ 2Tnf\y'^{\ + AM'^x'^/Q'^) 

+ 1^ 2 


in the i? ~ 0 limit. At high energies, the scalar and 
fermion evaluation of give the same numerical values 
for equal scalar and fermion masses because the energy 
loss is dominated by the small y region of integration. 

Our result for for f’s with a mass of 150 GeV is 
given by the upper curve in Fig. |21 For = 150 GeV, 
it is increasing with energy from ^ 2 x 10“® cm^/g at 
10® GeV to 3 X 10-® cm^/g at lO^^ GeV. 

The dashed curves in Fig. |21 represent the photonu¬ 
clear contribution to (3 for three stau masses: 50 GeV, 
150 GeV and 250 GeV, showing that f3 has a l/rri; depen¬ 
dence on the stau mass. This can be understood based 
on the argument which has been discussed in the con¬ 
text of electron-positron production in muon scattering 
by Tannenbaum in Ref. [9j. Qualitatively, 


F2 


+ ml 


/2 


( 12 ) 


where to first approximation, /2 is constant for small y. 
The parameter ml is 0.32 GeV^ for the Abramowitz et 
al. parameterization [^. We can approximate 


a 



Atto^ ^ 

g2(Q2+^2)/2 . 


(13) 


The integral over is approximately independent of 
^^‘1 depends on Qlain/^o- l’^ small y limit, 
~ mjy'^. The change of variables of the integral 
from y to z where 


z = 


mjy 

Too 


(14) 


yields a Jacobian factor of m^/mj multiplying the inte¬ 
gral. The stau mass dependence of /3 is given by this 
Jacobian factor. 


B. Stau Energy Loss due to Pair Production 

Next we consider e“'"e“ pair production in collisions of 
staus with nuclei. In this section, we show that the dom¬ 
inant contribution to /3 for the production of an e“'"e“ 
pair by a scalar charged particle interaction with a bare 
nucleus is the same, at high energies, as for the fermion 
case. We then use the results of Kokoulin and Petrukhin 
for fermion scattering M which include all of the dia¬ 
grams and incorporate a smooth transition between the 
unscreened to screened nucleus. 

This process is a very important contribution to the 
energy loss of muons and taus. In the case of staus, 
it is also important. There are several diagrams that 
contribute. They are shown in Fig.^ 

The diagrams can be classified into two categories. Di¬ 
agrams (al) and (a2) are “two photon” contributions 
to e+e“ pair production. Diagrams (bl), (b2) and (c) 
combine to make the gauge invariant set of “virtual 
bremsstrahlung” contributions to e“''e“ pair production. 
For taus instead of staus, the virtual bremsstrahlung 
graphs are of the type (bl) and (b2). The interference 
term between two photon (charge conjugation even) and 
virtual bremsstrahlung (G-odd) contributions vanishes in 
the integral. 

At relativistic energies, the two photon contribution 
dominates because of the presence of two quasi-real pho¬ 
tons. This is independent of the spin of the charged par¬ 
ticle (tau or stau) that emits one quasi-real photon. The 
two photon process leads to a factor proportional to the 
square of the classical electron radius Xe = e^/rUeC^, while 
the virtual bremsstrahlung contribution is proportional 
to the square of the classical heavy lepton radius [ia,a 
result that carries over to the heavy slepton case at high 
energy. The dominance of the two photon process was 
first demonstrated in Ref. m and has been studied in 
many contexts, includin g lep ton pair production in col¬ 
liding beam experiments , and more recently in heavy 
ion collisions [l4l| . 

Since the two photon diagrams dominate, we calculate 
only these two contributions. The cross section corre¬ 
sponding to these diagrams, in the notation of Kelnerp^. 
is given by: 


2„,4 


da = 


Z^a 


■[{Pi 




^ II 

kfQ4^ 
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One finds that G = G(a;+, a;_, 2/) is 


G 


ai 


bi 

Cl 


fj_ _ ^+^- \ _ L 

7TOg / ^ Q^x+x- 

x+x_ 

1 77j2^2 

{2xl + 2xl + l){D + XG) 

+2Z?^(x^ + x^ + 1) 

D(2a;^ + 2a;?. + 1) 

D(1+0 + AG 

+2a;+a;_[AG- (2 +A)£i] 


with the additional definitions 



—|■a;+a;_2A 
mi 


2 ( 1 - 2 /) ■ 


(18) 


(19) 

( 20 ) 


FIG. 4: Feynman diagrams for e''"e“ pair production by stau 
scattering from a nucleus with nuclear momentum transfer q. 
Diagrams (al) and (a2) dominate 

dp+dp_ dp2 dq ^ 

—0-0-0- b[q + ki - p+- p-} (15) 

p'± pI “0 

where is the velocity of the nucleus, and pi, p 2 , p+ 
and p- are the momentum for incoming stau, outgoing 
stau, (outgoing) position and electron, respectively. In 
addition, fci = pi— P 2 , Q = P+ + P- — ki, and = —q^. 
For scalar interactions, 

= (pi+P2)“(pi+P2)^ • (16) 

The quantity is the usual fermion trace involving the 
sum of the two diagrams with the two photon vertices. 
A series of separations and integrals can be performed so 
that the differential cross section can be written in the 
following form [T^: 

da = - y)dx+—dQ^ (17) 

o Trmg UJ 

'G(^x-i-^ X—^ Q , y) , 

where 


The quantities G and D characterize the incident 
particle as scalar or fermion: 

G = 0 D = 1 scalar case , 

G = 1 D = 1 fermion case . 

The resulting fd from diagrams (al) and (a2) of Fig.^ 
yielding Eq. IIHJ, is obtained for the scalar case and 
found to be the same as for the fermion case to within 
the numerical accurac y of our integration using the for¬ 
tran program VEGAS [I^ (here, a relative error of 10“"^). 
This can be understood explicitly from Eq. (18), where 
each factor of C, which differs for scalar and fermion 
cases, is multiplied by A. The quantity fd is dominated 
by small y, so A ^ 1. On a more qualitative level, one 
notes that the Gordon decomposition of the fermion cur¬ 
rent has a term proportional to {pi + P 2 )in the scalar 
coupling, which dominates in the high energy limit. 

Because the two photon ((al) and (a2) diagrams) pro¬ 
cess gives and at high energies, the difference be¬ 

tween fermion and scalar scattering to produce the e“'"e“ 
pair is negligible, in the results we show below we have 
used the complete fermion expression of Ref. M with 
the substitution ^ mj to describe for staus at 
high energy. 

The results for are shown by the labeled curve 

in Fig. [21 for mj = 150 GeV. The scaling with mass is 
shown in Fig.jSl where we show with solid lines for 


a;+ 


X- 
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FIG. 5: Stau range (solid) and time dilated decay length (dot- 
dashed line) for different values of '/F and for mj = 150 GeV 
and rrii = 250 GeV. The values of are chosen to make the 
lifetimes equal. The dashed line shows the parameterization 
of the range using Eq. (24-25). 


lO 



FIG. 6: Stau range (solid) and time dilated decay length (dot- 
dashed line) for = 150 GeV and m,;- = 250 GeV for x/F = 
10^ GeV. The dashed line shows the parameterization of the 
range using Eq. (24-26). 


mj = 50, 150 and 250 GeV. Arain, the 1/mj behavior 
of (3 appears. As noted in Ref. [3, the differential cross 
section (in y and Q^) can be approximated in the small y 
limit by a form similar to Eq. C3, where instead /2 —*■ 
<Je+e- 1 the photopair production cross section which is 
nearly constant, and toq —> 4me related to the threshold 
for electron-positron pair production. The same Jacobian 
factor, now ~ me/mj determines the mj dependence of 

opair 


C. stau Energy Loss due to Bremsstrahlung 


The bremsstrahlung contribution to /3 is much smaller 
than the pair production and photonuclear contributions 
by an additional factor of l/mp A textbook discussion 
of shows that 


da a 
dy y 


( 21 ) 


where d = lm~. The expression for /jt’rem 
has a dependence on the particle mass through ci, namely 
Ijm-. This is explicitly shown in Fig. |21for mj- = 150 
GeV, where the evaluation of fd follows from Ref. . 


III. STAU RANGE 


Our evaluation of the stau range follows the same pro¬ 
cedure for tau propagation outlined in Ref. 0. Briefly, 
we perform a one-dimensional Monte Garlo evaluation 
which treats the first two terms in the energy loss 


dE 

dX 


= a 


N 


-E 


f'Vc-U.t 

dyy 


N . 


dy A 


-E 


dyy 


' ycut 


da 
dy ’ 
( 22 ) 


continuously, and the third term stochastically. We 
use ycut = 10“^. We evaluate the survival probability 
P{E, Eq^ z) as a function of initial energy E, minimum 


final energy Eq and depth z. The average range is 


X{E,Eo) = J dzPiE,Eo,z) . 


(23) 


We note that the average range is not the same as the 
range of average energy loss, the solution to Eq. O US- 
Discussed below is a parameterization of an effective /3 
that, when substituted into the solution to Eq. o, gives 
the approximate average range. 

The average ranges of the f for two masses, mj = 150 
GeV and mj = 250 GeV are shown in Fig. El by the solid 
lines. For this figure and Fig. El the stau is propagated 
until its energy is reduced to Eq = 10^ GeV. Shown in 
Fig. El are the results for x/F = 10^ GeV for mj = 150 
GeV and vF = 2.15 x 10^ GeV for the larger mass. The 
two values of E are chosen to give the same lifetimes for 
mj- = 150 and 250 GeV. The time dilated decay length is 
long compared to the distance over which energy is lost 
by the stau. This is clear in that we can parameterize 
the range without reference to the decay length. A con¬ 
venient choice is to use the solution to Eq. o for the 
range of average loss with a constant /3, but to use in the 
solution an energy dependent (3. In terms of initial stau 
energy E and final stau energy Eg, we write 


X{E,Eq) 

= {a + PE)/{a + PEq) 

(24) 

p 

~ Pq + Pi\n{E/10^° GeV) 



+p 2 HEo/lO^ GeV) 

(25) 
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FIG. 7: Stau range (solid) for '/F = 10^ GeV and for mf = 
150 GeV and rrij- = 250 GeV, with Eq = 10® GeV. The dashed 
lines show the parameterization of the range following Eqs. 
(24-26) and the solid lines show the Monte Garlo results. 

with 

a = 2x 10“^ GeVcmVg 
/3 q = 5 X 10“® cm^/g(150 GeV/m-j-) 
f3i = 2.8 X 10“^° cm^/g(150 GeV/wj) 

/32 = 2 X 10“^° cm^/g(150 GeV/mj) . 

These parameters are obtained from a fit to the average 
range determined by the one-dimensional Monte Carlo 
program. In Fig. the dashed lines following the solid 
curves show the parameterized fit, Eqs. iI23J) and (ESI , 
with Eq = 10^ GeV. 

The average range of stau varies roughly linearly with 
stau energy E for lower energy E < 10® GeV in Fig. 
0 For nil = 150 GeV and Vf = lO’’ GeV, the range 
is about 40 km water equivalence (w.e.) for E = 10"* 
GeV, and about 3 x 10® km w.e. for E = 10® GeV. This 
behavior can be obtained from Eq. (1^ where for (3E <C 
a, X (X (E — Eq) / a. For higher energy, the average range 
varies logarithmically with energy: X oc \og{E/Eo)/P, 
roughly a few x 10"^ km w.e. 

The average range at low energy has little mass depen¬ 
dence since a is a nearly a constant as a function of mass. 
At higher energy, the larger the mass of the particle, the 
longer the average range due to the 1 /mj dependence of 

P- 

For y/F fixed at \/F = 10^ GeV, in the case of 
mj = 150 GeV, the low energy range is dominated by 
ionization energy loss, while for = 250 GeV, it is 
dominated by the decay length (see Fig.EI). For the lat¬ 
ter case, the average range at low energy still has the 
linear dependence on E since ■jct = {E/mi)cT. A mod¬ 
ification of Eq. ini) to account for the lifetime involves 
the substitution 

TTL~ 

a ^ ac = a -\ - - . (26) 

crp 



FIG. 8: Ratio of stau range using the parameterization of Eqs. 
(24-26) to stau range using a constant /3 = 0.8 x 10~^mTImj 
crrP/g for -\/F = 10^ GeV and mj- = 150 GeV (solid) and 
mj = 250 GeV (dashed line), and for -\/F = 2.15 x lO’’ GeV 
with mj- = 250 GeV (dot-dashed line), for Eq = 10® GeV. 

With Eq. ll^ and the energy dependent f3, the substi¬ 
tution oi a ac leads to a reasonable approximation to 
the range, as seen in Fig. El Numerically, for the three 
range curves shown in Figs. El and El the parameterization 
agrees to within ~ 20% over the whole energy range. In 
fact, for the curves in Fig. El the ratio of parametrization 
of the range to the Monte Carlo evaluation of the range 

deviates from 1 only by —8% -1-1% from 10"^ — 10^^ 

GeV initial stau energy. For the shorter lifetime shown 
in Fig. El the ratio of the parameterization to the Monte 
Carlo range deviates from 1 by —5% —1-15% except for 
10"^ — 10® GeV, where the deviation reaches as much as 
- 22 %. 

At higher energy thresholds Eq, the effective /3 in¬ 
creases, requiring the P 2 term in Eq. (123) • In Fig. 13 
we show the stau range and its parameterization for 
Eq = 10® GeV. The average range is shorter than the 
case of Eq = 10® GeV due to the higher energy thresh¬ 
old. The ratio of the parameterization of the range to 
the Monte Carlo evaluation goes between 1.03 — 0.88 for 
mj- = 150 GeV, and between 1.11 — 0.93 for mj- = 250 
GeV. 


IV. DISCUSSION 

Albuquerque, Chacko and Burdman proposed to look 
for evidence of supersymmetry in a class of models 
by looking for pair-produced staus with nearly parallel 
tracks Q. They concluded that there is a potential sen¬ 
sitivity of neutrino telescopes to a range of parameter 
space in certain low energy SUSY breaking models. Their 
evaluation of event rates of upward-going staus depends 
in part on the effective area of the detector and the range 
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of the staus. 

To first approximation, the energy loss of the stau can 
be characterized by constant a and This is what 

was done in Ref. Q, where Pr — 0.8 x 10“® cuP/g for 
the tau is rescaled to Pf^c — 9.5 x 10“® (150 GeV/wj-) 
cm^/g according to the factor of ijit/itij. With our new 
estimate of the stau range, parameterized by an energy 
dependent /3, the effective range is increased relative to 
the range calculated using Pf c- 

Fig.ia shows the ratio of the range using the energy 
dependent P compared to the range calculated using 
the constant Pf^c, scaled appropriately for the different 
masses, for initial stau energies between — 10^^ GeV 
and a minimum stau energy of Eg = 10^ GeV. The ratio 
can be as much as a factor of two for initial stau energies 
in the range of 10^ — 10® GeV. In this energy range, the 
P term in Eq. CJ dominates over the ionization term (a) 
and the lifetime, but the effective P is not as large as the 
rescaled constant /3f in this energy range. The enhance¬ 
ment of the range by a factor of two is significant because 
it substantially improves the potential of neutrino tele¬ 
scopes for detecting evidence of stau production. 

For Eq = 10® GeV, a similar comparison between the 
Monte Carlo range and the range using is obtained. 
The ratio of the range determined by Monte Carlo prop¬ 
agation to the range evaluated with Eq. o using Pf,c is 
a factor of ^ 2 for E = 10^ — 10® GeV for mf = 150 GeV 
and a factor of about 1.6 for the same energy range with 
mi = 250 GeV for y/F = 10^ GeV. 

Estimates which scale the energy loss P by 1/m? 
dramatically overestimate the stau range at high ener¬ 


gies. Using the P of Ref. |^, /3 = 3.9 x 10“®(m^/mj')^ 
cm^/g, one finds that for an initial stau energy of 10® 
GeV, one has overestimated the range by about 40%, 
and at 10^ GeV by a factor of 4.8. The overestimate 
rapidly increases with energy. It is a factor of ^ 30 for 
an initial energy of 10® GeV, a factor of ^ 1600 for an 
initial energy of 10^^ GeV. 

In summary, for the dominant photonuclear and e“'"e“ 
pair production processes, one expects qualitatively that 
the /3 of a charged scalar particle (stau) scales as the in¬ 
verse mass of the scalar particle. We have shown this 
result quantitatively by explicit calculation. Incorporat¬ 
ing the electromagnetic energy loss of the stau and the 
long lifetime following from the supersymmetry breaking 
scale, we used a Monte Carlo to evaluate the range. The 
range determined by Monte Carlo evaluation is longer 
than the range of average energy loss using Eq. o and 
a constant /3f,c- Our energy dependent parametrization 
of an effective P in Eq. (1^ for the stau range reasonably 
approximates the Monte Carlo results. Our evaluation of 
the stau range suggests that previous estimates Q of up¬ 
ward stau pair event rates, induced by supersymmetric 
particle production by upward-going neutrinos, may be 
enhanced by a factor of ~ 2 in some energy ranges. 
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